Violations of local realism by two entangled quiVits are stronger than for two qubits 
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Tests of local realism vs quantum mechanics based on Bell's inequality employ two entangled 
qubits. We investigate the general case of two entangled quiVits, i.e. quantum systems defined in an 
iV-dimensional Hilbert space. Via a numerical linear optimization method we show that violations 
of local realism are stronger for two maximally entangled quiVits (3 < N < 9), than for two qubits 
and that they increase with N. The two quiVit measurements can be experimentally realized using 
entangled photons and unbiased multiport beamsplitters. 
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John Bell M has shown that no local realistic mod- 
els can agree with all quantum mechanical predictions 
for the maximally entangled states of two two-state sys- 
tems (qubits). After some years researchers started to 
ask questions about the Bell theorem for more compli- 
cated systems. The most surprising answer came from 
the GHZ theorem || : for three or more qubits the conflict 
between local realism and quantum mechanics is much 
sharper than for two qubits. The other possible extension 
are entangled states of pairs of A^-state systems, quTVits, 
with N > 3. First results, in 1980-82, suggested that the 
conflict between local realism and quantum mechanics 
diminishes with growing N |j| . This was felt to be in con- 
currence with the old quantum wisdom of higher quan- 
tum numbers leading to a quasi-classical behavior. How- 
ever, that early research was confined to Stern-Gerlach 
type measurements performed on pairs of entangled 
spins ||. Since operation of a Stern-Gerlach device de- 
pends solely on the orientation of the quantization axis, 
i.e. on only two parameters, devices of this kind cannot 
make projections into arbitrary states of the subsystems. 
That is, they cannot make full use of the richness of the 
iV-dimensional Hilbert space. 

In early 1990's Peres and Gisin 0] considered certain 
dichotomic observables applied to maximally entangled 
pairs of quA^its. They showed that the violation of lo- 
cal realism, or more precisely of the CHSH inequalities, 
survives the limit of N — > oo, but never exceeds the vi- 
olation by two qubits, in agreement with Cirel'son limit 
|1, i.e. it is limited by the factor of \/2. Therefore, the 
question whether the violation of local realism increases 
or not with growing TV for general observables was still 
left open. 

To answer this question it is necessary first to adopt 
an objective measure of the magnitude of violation of 
local realism. To this end, consider two quAlt systems 
described by mixed states in the form of 

Pn(F n ) = F NPnmse + (1 - F N )\H>» ax )(*Z ax \, (1) 

where the positive parameter Fn < 1 determines the 
"noise fraction" within the full state, p no ise = -ffili an d 
Y^max) i s a maximally entangled two quTVit state, say 

1 N 

K ox } = -t=y, i™umb- (2) 

m— 1 

In (||) \iti)a (|to)b) describes particle A (B) in its mode 
to. One has x {m\Tn>')x = o~m,m', with x — A,B. The 
threshold maximal F N lax , for which the state pn(Fn) 
still does not allow a local realistic model, will be our 
value of the strength of violation of local realism. The 
higher F^ ax the higher noise admixture will be required 
to hide the non-classicality of the quantum prediction. In 
experiments the visibility parameter V, effectively equiv- 
alent to 1 — Fn, is the usual measure of the reduction of 
interferometric contrast (visibility). 



We shall study the case of two observers Alice and 
Bob performing measurements of local non-degenerate 
observables, each on her/his quAIt of an entangled pair in 
the state pn{Fn)- Let us imagine that Alice can choose 
between two non-degenerate observables A\ and Ai, and 
that each observable is defined such that it has the full 
spectrum characterized by all integers from k = 1 to N. 
Bob can choose between B\ and B 2 , both with the same 
spectrum as above (/ = 1,2, ...,N). Thus, the observers 
can perform 2x2 mutually exclusive global experiments. 
The quantum probability for the specific pair of results, 
k for Alice and I for Bob, provided a specific pair of lo- 
cal observables is chosen, Ai by Alice and Bj by Bob, 
will be denoted by Pp (k; l\Ai, Bj). Quantum mechan- 
ics makes predictions for the complete set of AN 2 such 
probabilities, and nothing more. 

The hypothesis of local hidden variables tries to go 
beyond. The basic assumption there is that each parti- 
cle carries a probabilistic or deterministic set of instruc- 
tions how to respond to all possible local measurements it 
might be subject to. Therefore local realism assumes the 
existence of non-negative joint probabilities involving all 
possible observations from which it should be possible to 
obtain all the quantum predictions as marginals (see, e.g. 
||, 0). Let us denote these hypothetical probabilities by 
P HV (k, to; I, n\Ai, A2, B\, B2), where k and to, represent 
the outcome values for Alice's observables (I and n for 
Bob's). In quantum mechanics one cannot even define 
such objects, since they involve mutually incompatible 
measurements. The local hidden variable probabilities 
for the experimentally observed events, k (to) by Alice 
measuring Ai (A2), and I (n) by Bob measuring B\ (B2), 
are the marginals 

P HV {k; l\A u BO = £ m £„ P HV {k, to; I, n), 
P HV (k; n\A u B 2 ) = £ m £ ; P HV (k, to; I, n), 
P HV (m; l\A 2 , B{) = £ fc £„ P HV (k, to; I, n), 
P HV (to; n\A 2 , B 2 ) = Ei P HV (*, ™; h »), 

(3) 

where P HV (k,m;l,n) is a short hand notation for 
P HV (k,m;l,n\A 1 ,A 2 ,B 1 ,B 2 ). The AN 2 equations (§) 
form the full set of necessary and sufficient conditions 
for the existence of local realistic description of the 
experiment, i.e., for the joint probability distribution 
P HV (k,m;l,n). The Bell Theorem says that certain pre- 
dictions by quantum mechanics are in conflict with the 
local hidden variable model (|[). Evidently, the conflict 
disappears when enough noise is added, as in the state 
(|l|), since that noise has a local realistic model. There- 
fore a threshold F N lax exists below which one cannot 
have any local realistic model with P HV (k; l\Ai, Bj) = 
Pp^ 1 (k; l\Ai, Bj). Our goal is to find observables for the 
two quA^its returning the highest possible critical F N lax . 
Up to date, no one has derived Bell-type inequalities 
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that are necessary and sufficient conditions for (p3) to 
hold, with the exception of the A — 2 case (see 0). 
However there are numerical tools, in the form of the 
very well developed theory and methods of linear opti- 
mization, which are perfectly suited for tackling exactly 
such problems J8|. 

The quantum probabilities, when the state is given by 
(|l|), have the following structure 

- &F N + (1 - F N )pQ M {k; l\Ai, Bj), (4) 

where P® M (k; l\Ai, Bj) is the probability for the given 
pair of events for the pure maximally entangled state. 
The set of conditions (||) with Pp M (k; l\Ai, Bj) replac- 
ing P HV (k; l\Ai, Bj) imposes linear constraints on the iV 4 
"hidden probabilities" P HV (k, m; I, n) and on the param- 
eter Fn, which are the nonnegative unknowns. We have 
more unknowns (A 4 + 1) than equations (4 A 2 + 1, with 
the normalization condition for the hidden probabilities), 
and we want to find the minimal Fn for which the set 
of constraints can still be satisfied. This is a typical lin- 
ear optimization problem for which lots of excellent algo- 
rithms exist. We have used the state-of-the-art algorithm 
HOPDM 2.30. (Higher Order Primal Dual Method) @. 
It is important to stress that for cross-checking four in- 
dependently written codes were used, one of them em- 
ploying a different linear optimization procedure (from 
the NAG Library). 

We were interested in finding such observables for 
which the threshold F/v acquires the highest possible 
value. To find optimal sets of observables we have used 
a numerical procedure based on the downhill simplex 
method (so called amoeba) [jl0| . If the dimension of the 
domain of a function is D (in our case D = An, where n 
is the number of parameters specifying the nondegener- 
ate local observables belonging to a chosen family), the 
procedure first randomly generates D + 1 points. In this 
way it creates the vertices of a starting simplex. Next it 
calculates the value of the function at the vertices and 
starts exploring the space by stretching and contract- 
ing the simplex. In every step, when it finds vertices 
where the value of the function is higher than in others, 
it "goes" in this direction (see e.g. |1C(]). 

Let us now move to the question of finding a family of 
observables, which returns critical -Fjv's that are above 
the well known threshold for the two qubit case, 1 — 
As it was said earlier, and was confirmed by our 
numerical results, Stern-Gerlach type measurements are 
not suitable. More exotic observables are needed. 

First we discuss how experiments on two entangled 
quAits might be performed. In view of the unavail- 
ability of higher spin entanglement it is fortunate that 
quAit entanglement can be studied exploiting momen- 
tum conservation in the many processes of two-particle 



generation, most notably in the parametric down conver- 
sion generation of entangled photon pairs. This results 
in strong correlations between the propagation directions 
of the particles in a pair. One can then submit A spatial 
modes of each particle to a multiport beamsplitter 

Application of multiports in the context of quantum 
entanglement has been first discussed by Klyshko Jll| . 
Proposals of Bell experiments with the multiports were 
presented in |Q, and further developed in iQ. Multi- 
port devices can reproduce all finite dimensional unitary 
transformations for single-photon states Q], therefore 
they are characterized by A 2 — 1 real parameters. 

In order to limit computer time we restricted our anal- 
ysis to unbiased multiports more specifically to 
Bell multiports. Unbiased multiports have the property 
that a photon entering into any single input port (out 
of the A), has equal chances to exit from any output 
port. In addition, for Bell multiports the elements 
of their unitary transition matrix, XJ N , are solely pow- 
ers of the N-th root of unity 7 at = exp(i2ir/N), namely 

Let us now imagine two spatially separated experi- 
menters who perform the experiment of FIG. 1. (de- 
scribed in the caption). The initial maximally entangled 
state (^|) of the two quAits can be prepared with the aid 
of parametric down conversion (see (l| ) . The two sets of 
phase shifters at the inputs of the multiports (one phase 
shifter in each beam) introduce phase factor e 1 ^- 4 " 1 "^) 
in front of the m-th component of the state (||), where 
4>™ and denote the local phase shifts. 




FIG. 1. The experiment of Alice and Bob with entangled 
quTVits. Each of their measuring apparata consist of a set of N 
phase shifters just in front of an 2N port Bell multiport, and 
N photon detectors Dk,Di (perfect, in the gedanken situation 
described here) which register photons in the output ports of 
the device. The phase shifters serve the role of the devices 
which set the free macroscopic, classical parameters that can 
be controlled by the experimenters. The source produces a 
beam-entangled two particle state. 

Each set of local phase shifts constitutes the interfer- 
ometric realizations of the "knobs" at the disposal of 
the observer controlling the local measuring apparatus, 
which incorporates also the Bell multiport and N detec- 
tors. In this way the local observable is defined. Its 
eigenvalues refer simply to registration at one of the A 
detectors behind the multiport. The quantum prediction 
for the joint probability Pp M (k, I) to detect a photon at 
the fc-th output of the multiport A and another one at 
the l-th output of the multiport B is given by (l| : 
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(^) (N + 2(1 - F N ) YZ>n «* (*33 - *£)) > ( 5 ) 



where = <j>J + $g + [m(k + 1 - 2)] ^ . The counts at a 
single detector, of course are constant, and do not depend 



upon the local phase settings: Pp (k) = Pp (I) = 
1/N. 

The numerical values of the threshold Fn are given in 
fig. 2. 



0.31879 



0.31793 
0.31674 

0.31511 I 

0.31285 f 



0.30945 



0.30385 



0.29289 



0.320 
0.316 
0.312 
0.308 
0.304 
0.300 
0.296 
0.292 
0.288 



N 

of pure noise admixture to 
a maximally entangled two quiVit system, such that a local 
realistic explanation still cannot be upheld. For smaller noise 
fractions a conflict arises between quantum mechanics and 
local realism. The result for N — 2 agrees with the standard 
threshold of 1 \=. 

v 2 



It is evident, that indeed two entangled quTVits vio- 
late local realism stronger than two entangled qubits, and 
that the violation increases monotonically with TV. It is 
tempting to contemplate the limit of N — > oo. While ob- 
viously the values of FJ^ ax seem to saturate, at present 
we cannot give a definite asymptotic value. 

A few words of comment are needed. One may argue 
that because of the rather large number of local macro- 
scopic parameters (the phases) defining the function to 
be maximized with the amoeba we could have missed 
the global minimum. While this argument cannot be 
ruled out in principle, we stress that in that case the ul- 
timate violation would even be larger. This would only 
strengthen our conclusion that two entangled quiVits are 
in stronger conflict with local realism than two entangled 
qubits. 

Based on the numerical results, i.e. the values of the 
optimal phase settings, and on the structure of the local 



hidden variable model for F™ ax , an algebraic calcula- 
tion was performed Jl5| showing that for the two qutrits 
(N = 3) experiment the exact value for F™ ax is n ~~ • 
One should also mention that for two spin 1 particles in a 
singlet state observed by two Stern Gerlach apparatuses 
our method gives F^ G = 0.1945, which is much smaller 
than 1 ' con fi rmm g that such measurements are not 
optimal in the sense of leading to maximal possible vio- 
lations of local realism. 

An important question is whether unbiased Bell multi- 
ports provide us with a family of observables in maximal 
conflict with local realism. For a check of this question we 
have also calculated the threshold value of F3 for the case 
where both observers apply to the incoming qutrit the 
most general unitary transformation belonging to a full 
SU(3) group (i.e. we have any trichotomic observables on 
each side). Again we have assumed that each observer 
chooses between two sets of local settings. However, in 
this case each set consists of 8 local settings rather than 
the three (effectively two) in the fritter case. The re- 
sult appears to be the same as for two fritters. While 
this might suggest that for N — 3 Bell multiports are 
optimal devices to test quantum mechanics against local 
realism, this needs to be further investigated. 

It is interesting to compare our results with the limit 
for the non-separability of the density matrices (|l|). The 
critical minimal Fn for which a density matrix ([!]) is 
separable is jM^ (see ]l6|). The fact that this limit is 
always higher than ours indicates that the requirement of 
having a local quantum description of the two subsystems 
is a much more stringent condition than the requirement 
of admitting any possible local realistic model. 

It will be interesting to consider within our approach 
different families of states, generalizations to more than 
two particles, extensions of the families of observables 
and to see if more than two (e.g. Ai, A2, A3) experi- 
ments performed on either side can lead to even stronger 
violations of local realism. The questions concerning the 
critical Fn are also important in the attempts to gener- 
alize Ekert's quantum cryptographic protocol to qutrits 
and higher systems [ p7[ . 

Our method is numerical, and is based on linear opti- 
mization. It is a development of the approach of ||. The 
exploding (with N) difficulty of approaching this type of 
problems via algebraic-analytical methods (generalized 
Bell inequalities, via the Farkas lemma, etc.) has been 
exposed in 0. 

It will certainly be fascinating to see laboratory real- 
izations of the experimental schemes discussed here. 
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